Analytical solutions for two heteronuclear atoms in a ring trap 
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We consider two heteronuclear atoms interacting with a short-range 8 potential and confined in 
a ring trap. By taking the Bethe-ansatz-type wavefunction and considering the periodic boundary 
condition properly, we derive analytical solutions for the heteronuclear system. The eigen-energies 
represented in terms of quasi-momentums can then be determined by solving a set of coupled 
equations. We present a number of results, which display different features from the case of identical 
atoms. Our result can be reduced to the well-known Lieb-Liniger solution when two interacting 
atoms have the same masses. 
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I. INTRODUCTION 

Motivated by recent experimental studies on heteronu- 
clear atoms, mixtures of atomic gas with different masses 
have recently attracted lots of attention. Stable mixtures 
of heteronuclear atoms, for examples, 40 K and 6 Li, 40 K 
and 87 Rb, and 6 Li and 23 Na, have already been real- 
ized in experiments [H-Q- The heteronuclear quantum 
gas mixture offers a wide range of possibilities for explor- 
ing novel phases and new quantum effects. Significant 
efforts have been made on research topics such as su- 
perfluid properties and pairing mechanisms of heteronu- 
clear fermions 043 an d heteronuclear molecule forma- 
tion 0, E3-E3- 

The mass ratio between heteronuclear atoms adds a 
new degree of freedom for exploring physical effects as- 
sociated with the mass imbalance. In general, heteronu- 
clear atomic systems are expected to behave quite dif- 
ferently than the equal-mass counterpart ED, [Hj]. 
To gain deep insight into properties of heteronuclear sys- 
tems, some analytical solutions are especially important. 
While most of theoretical works on the many-body het- 
eronuclear systems are within the mean-field theory, few- 
body systems can be treated much more precisely within 
nonperturbative microscopic frameworks 14]. Particu- 
larly, the two-body problem is vital to understand the ba- 
sic physical phenomenon from the beginning of quantum 
mechanics. The two-atom system composed of identical 
atoms is exactly solved as one can separate the center- 
of -mass and relative motions and thus the two-particle 
system turns to a one-body problem p^| - [T9j . The prob- 
lem with different atoms becomes complicated and gen- 
erally is hard to get analytical solutions as it does not 
separate in center-of-mass and relative coordinates [l2j]. 
Recently, the scattering of heteronuclear atom pairs in a 
one-dimensional optical lattice has been studied by using 
the method of Green's function 20]. 

In this article, we study the problem of two unequal- 
mass atoms with (^-potential interaction in a one- 



dimensional ring trap with periodic geometry. The pe- 
riodic geometry has been experimentally realized as a 
circular wave guide or toroidal trap [2l|, |22j . For the two- 
atom system with (^-potential interaction, it is convenient 
to treat the interaction term by separating center-of-mass 
and relative motion. However, it is not convenient to ap- 
ply the periodic boundary condition by using separated 
center of mass and relative coordinates. When the peri- 
odic boundary condition is taken into account, the rel- 
ative motion is actually coupled to the center-of -mass 
motion. In the present work, we analytically solve the 
problem of interacting two heteronuclear atoms by taking 
the Bethe-ansatz-type wavefunction and considering the 
periodic boundary condition properly. While the momen- 
tum of center-of-mass motion is a good quantum number, 
the relative momentum is coupled to the center-of-mass 
momentum due to the periodic boundary condition. Our 
results are consistent with the exact Bethe-ansatz solu- 
tions of the Lieb-Linger model when two atoms have the 
same mass 18]. With the analytical solution, we obtain 
the energy spectrum, density distribution and momen- 
tum distribution for different mass ratio. Our results in- 
dicate that they display different features from the iden- 
tical two-atom system. 

The paper is organized as follows: In Sec. II, we intro- 
duce the model and present a detailed derivation of the 
analytical solution of the interacting heteronuclear sys- 
tem. In the limit of equal mass, our analytical result 
is compared to the Bethe-ansatz solution of two-particle 
Lieb-Liniger model. In Sec. Ill, we present the energy 
spectrum of heteronuclear atoms. In Sec. IV, we show 
the density distribution of relative coordinate and the 
momentum distribution of heteronuclear atoms. A sum- 
mary is given in the last section. 



II. THE MODEL AND ITS SOLUTIONS 

We consider an interacting two-atom system in a ID 
ring trap described by the Hamiltonian 
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H=- 
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2mi dx\ 2to2 dx 2 . 



gS{xi-x 2 ), (1) 
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where g is the interaction strength and the general case 
with different masses mj and m 2 is considered. The ge- 
ometry of the ring trap enforces the periodic boundary 
condition 

9(x 1 +L,x 2 ) = V(xi,x 2 + L) = *(si,x 2 ) (2) 

to the eigenfunction ty(x\,x 2 ) of Hamiltonian H. In- 
troducing relative and center-of-mass coordinates x — 
x\ — X2, X — {m\Xi + 1TI2X2) / (nil +"12), or alternatively 
x\ = X H ^r 2 — x, xo = X — x, we can rewrite 

1 mi+rri2 ' z mi+m2 ' 

the Hamiltonian Eq.(JTJ as the separated form 
H = H x (X) + H x (x) 



h 2 d 2 h 2 d 2 



2MdX 2 2u,dx 2 ' *" y ~ n ™ 

where fi — mim 2 /(mi + 771.2) and M = mi + m 2 . The 
wavefunctions of the Hamiltonian can take the form of 
iff(X,x) — e lKX ip(x), where e lKX is the eigenstate of 
Hx{X) and the eigenstate of H x (x) has the form of 

ip(x) = (A+e tkx + B+e- tkx )6(x) 

+(A-e ikx + B_e- lkx )9(~x), (4) 

where A± and B± are coefficients to be determined. It 
is straightforward to get the eigenenergy given by 



E = —K z 
2M 



2[i 



-k 



(5) 



Although formally the center-of-mass and relative parts 
are separated, we shall show that the relative quasi- 
momentum k is actually coupled with the total momen- 
tum K due to the periodic boundary condition ([2]). 

The continuous condition of the wavefunction Q gives 
the following restriction to the coefficients: 



A++B+ = A. 



B- 



(6) 



By integrating the eigen-equation H'^>(X 1 x) = E^(X, x) 
from the negative infinitesimal — e to the positive in- 
finitesimal e, we get 



d 



2/i \dx 
which leads to 



d_ 

dx^ 



¥>(*),=_« + MO) =o, 



ik (A + -B+- A_+B_) = J^g (A+ + B+ + A-+B-). 

(7) 

The Eq.© and Eq.© can not uniquely determine all 
the parameters. Additional relations can be obtained 
by considering the periodic boundary condition (|2"|). In 
order to apply the periodic boundary condition, it is more 
convenient to work in the coordinate space of x\ and x 2 . 
In terms of xi and X2, the total wavefunction *£(xi, x 2 ) = 



,JKX 



ip(x) can be also represented as 



^ A+e iktxt+ik 2 x 2 + B + e< Xl+lk ^ X2 )e(xi - x 2 ) + 
( A _ e ik lXl +ik 2X2 + B_e ik 'i Xl+ik 2 X2 )6(x 2 - xi), (8) 



where k\ = Kmi/(mi +m 2 ) + k, k 2 = Km 2 / (mi+m 2 ) - 
k, k[ — Kmi/(mi+m 2 ) — k, and k' 2 = Km 2 / (m\-\-m 2 )- 
k. The periodic boundary condition then enforces 



A + e lklL 


= A 




= B 


A^e lk2L 


= A 


B_ e ik 'i L 


= B 
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(9) 
(10) 
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By multiplying Eqs. © and ((TTJ) or Eqs (TTOjand (fHI) . 

we get e lKL = 1, which implies 



KL = 2im 



(13) 



with n the integer. Combining Eqs. © (ITU)) and ©, we 
can represent B + , A_, £?_ in terms of A + . Substituting 
them into Eq. ((7|), we can derive the following relation 



k— = 

cos 



sin (kL) 



\ m 1 +m 2 



cos (kL) 



(14) 



From Eq. ([T4l) . it is clear that k is coupled with the K 
whereas K is uniquely determined by Eq. (ll3[) . By solving 
Eqs. (IT3|) and ([14]) . we can get solutions of k and K and 
therefore eigenvalues of the system by Eq. ([5]) . 

The above equations can be also represented in terms 
of the quasi-momentum ki and k 2 via K = ki + k 2 and 
k = (iri2ki — m\k 2 )/(m\ + 1112)- In terms of k\ and k 2 , 
the eigenenergy is given by 



E : 



2m\ 



2m 2 



(15) 



In order to compare to the case with the same masses, 



we introduce g = ^-c and rewrite Eq. ([14")) as 



1\ k 
1 + - - 

ale 



sin (kL) 



cos [wk) -cos(fcL) 



(16) 



where a = m 2 /mi is the mass ratio. When the mass ratio 
a = 1, our results (fP3")) and (fi"4]) should be consistent with 
the well-known Bethe-ansatz solution of the Lieb-Liniger 
model [l8j ]. Using the above relation (fl~6| and taking 
a = 1, we can get 



2k- 



ic 



2k 



e ikL _ CQS {KL/2) 
cos (KL/2)-e~ ikL ' 



Since cos(KL/2) = cos(7rn) = ±1, the above equation 
reduces to 



h - k 2 



ki 



= ±e 



IC 



where + (or — ) corresponds to n =even (or odd). By 
using the relation e lkL — e < k i- k ?) L / 2 — e ikiL e -iKL/2 _ 
±e lfcli , we then have 



a ik\L 



fcl 



IC 



ki — k 2 — ic ' 



(17) 
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which is just the Bethe-ansatz equation of the two- 
particle Lieb-Liniger model 18]. Another Bethe-ansatz 
equation can be obtained via e lk2L = e~ %klL . When 
mi = we have k[ — ki and k' 2 = ki. From Eqs. 
(fT0|) and ([6]), we can also get v4_ = B + and B_ = A + . 
Therefore the wavefunction also has the same form of 
the Lieb-Liniger model [l8[ , which is invariant under the 
the exchange of x\ and X2 ■ We note that, for the gen- 
eral case with mass imbalance, the wavefunction is not 
invariant under the exchange of x\ and X2 since k[ =/= k% 
and kl 2 ^k\. 



III. THE ENERGY SPECTRUM 
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Without loss of generality, only the case of a > 1 shall 
be considered. For convenience, in the following calcu- 
lation we shall take h = 1 and 2mi = 1. In the units 
of energy h 2 j(2m\) = 1, the energy expression Eq.([5]) or 
Eq. ([T5l) becomes 



E = (1 + a^K 2 + (1 + l/a) k 2 



(18) 



or E = k\ + (l/a)k 2 . Since the total momentum K 



is a conserved quantity determined by Eq. (|13|) . the en- 
ergy spectrum of the system E can be characterized by 
different K. The ground state corresponds to K = 0. 
Given K — 2im/L, k can be determined by solving 
Eq. (fT4| or Eq. (fTBj) . For the repulsive case with c > 0, 
Eq. (TT51) has only real solutions which describe the scat- 
tering states. However, for the attractive case with c < 0, 
Eq. (|16|) may have imaginary solutions corresponding to 
the bound states of the attractive system. 

First we shall consider the case with the total momen- 
tum K — 0. For K = 0, Eq. dTBl) becomes 



1 \ k 
1 + - - = 



sin (kL) 
a J c 1 — cos (kL) 

With similar steps as the equal-mass case, we can get 



ikL 



IC 



(19) 



with c' = c/(l + i). For c > 0, k has only real solutions. 
It is convenient to solve the real solution of k from the 
logarithm form of Eq.(fT9|) given by 



(20) 



kL = 2w(I - 1/2) - 2arctan(A;/c'), 



where / takes integer. By numerically solving the k so- 
lution, we can get the energy E = (l + —J k 2 for states 
with K = 0. The ground state energy corresponds to the 
solution of Eq.([20]) with 1 = 1. For c < 0, Eq.(jl9]) has 
imaginary solution k = iX. The real solution A can be 
obtained by solving 



-XL 



A 



FIG. 1: (Color online) The dimensionless energy spectrum 
versus 7 for K = 0, a = 2.175 (dashed line) and a = 1 (solid 
line). For 7 < 0, there exists a bound state. 



(21) 



which is obtained by inserting k = i\ into Eq. ([T9|) . From 
Eq. ([2~T|) . one can also check that A has only solution for 
d < 0. The energy of the bound state can be obtained 
via£ = -(1 + 1/a) A 2 . 

To give a concrete example, we consider a two-atom 
system composed of 40 K and 87 Rb with the mass ration 
a = 2.175. In Fig.l, we plot the energy spectrum for 
states with K = 0. The energy spectrum for a = 2.175 
is represented by dashed lines. The energy spectrum of 
the equal-mass system is also given for comparison. Here 
we use the dimensionless energy e( 7 ) = E/(Np 2 ), where 
p = N/L and N = 2 is the atom number. It is shown 
that the ground energy for the system of a > 1 has lower 
energy than that of the system with a = 1. Actually, all 
the other spectrum corresponding to quantum number / 
for a > 1 is below its correspondence for a = 1. 

Now we consider states with the nonzero total momen- 
tum, for example, K = ±^-, which correspond to excited 
states of the system. When K ^ 0, Eq. (fT6]) can not be 
simplified to a similar form as Eq. (TT9")) or Eq. ([2T)]) for an 
arbitrary a. Nevertheless, formally we have 

, k r kL o KL 

kL = 2tt1 — 2 arctan 2 arctan tan — cot — ; , 

d L 2 2(l + a) J ' 

(22) 

where / takes integer. When K = 0, it reduces to 
Eq. (120[) . The relative momentum can be obtained by 
numerically solving Eq- lfTB]) or (|22"T) with K — ± 2 f- and 

the energy is given by E — (1 + a) -1 (x") 2 + (l + ^) k 2 . 
The energies for K = 2tt/L and K = —2tt/L are degen- 
erate. In Fig. 2, we display energy spectrums of states 
with K = ±2ir/L for a = 2.175 and a = 1. It is shown 
that the energy spectrum for a > 1 is always below the 
corresponding spectrum for a = 1. 

The effect of mass ratio on energy spectrum can be un- 
derstood from two aspects. First, we can see that coeffi- 
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FIG. 2: (Color online) The dimensionless energy spectrum 
versus 7 for K = ±2n/L, a = 2.175 (dashed line) and a — 1 
(solid line). 



cients in Eq.QlSp decrease as the mass ration a increases. 
On the other hand, while K in Eq. tfLSl) is independent 
of the mass ratio, k is related with a via the nonlinear 
equation ([To]) . It is not easy to analytically analyze the 
influence of the mass ratio on k from Eq. (|16p . However, 
for the case of K = 0, from Eq. (fl9)) or Eq. (j20|) . we notice 
that the mass ratio a affects the solution of k by effec- 
tively changing the interaction strength via d and thus 
k increases with the increase in a. 

For the attractive case with c < 0, it is worthy of in- 
dicating that there exists a bound state for each set of 
states characterized by K when the interaction exceeds a 
threshold c t . The threshold of j t = c t/p for the appear- 
ance of bound state is given by 



1 I 1 

It = r 1 + - 

2 \ a 



cos 



KL 
T+a 



- 1 



(23) 



For K = 0, we have 7$ = 0. This implies that the atoms 
form bound state for arbitrary weak attraction 7 < 0. 
For K = 2tt/L, we have 74 < 0. In this case, atoms can 
form bound state only for 7 < 7t- 



IV. DENSITY AND MOMENTUM 
DISTRIBUTION 

To show the influence of the interaction on the wave- 
function, we plot the ground state density distribution 
for the relative motion defined as p(x) — ip*(x)tp(x) in 
Fig. 4. The density distribution displays a cusp at x = 
due to the repulsive interaction when two atoms coin- 
cide. The relative motion density distribution at x = 
decreases as the interaction strength increases. The den- 
sity distribution at x = for the case of a > 1 is below 
that of the equal-mass case as a consequence of increasing 



FIG. 3: (Color online) The probability density distribution of 
relative motion for the ground state with 7=1 and a = 1 
(dash dot line), 7 = 1 and a = 2.175 (dashed line), 7 = 5 
and a = 1 (short dashed line), 7 = 5 and a = 2.175 (dashed 
dotted dot line), 7 = 500 and a — 1 (short dashed dotted 
line), 7 = 500 and a = 2.175 (red dotted). 



c! effectively. In the limit of c — > 00, the wave function 
Sfr(xi,X2) vanishes whenever x\ — Xi. 

The momentum distribution for atom with mass mi 
can be represented as 



where ^>(pi,p2) is the Fourier transformation of 
^(xi, x<i). The definition of momentum distribution 
n 2 (p) for the atom with mass m 2 is similar. For the equal- 
mass case, the momentum distribution for atom 1 or 2 
should be the same since the wave function ^{xx^x?) is 
invariant under the exchange of X\ and x 2 - However, the 
unequal-mass atomic system will show different features 
since the wave functions are not necessary to be sym- 
metrical for exchanging two unequal mass atoms except 
of the case of K — 0. For K = 0, we find that the wave- 
function is invariant under exchange of X\ and x 2 even for 
the unequal-mass case, and thus we have ni(p) — 7ia(p)- 
In Fig. 4, we display the ground state momentum distri- 
bution for the unequal-mass system with a = 2.175 and 
the equal-mass system. An obvious peak around the zero 
momentum position can be observed. With the increase 
in the interaction strength, the momentum distribution 
tends to spread out widely and the height of the peak 
decreases. There is no obvious difference for momentum 
distributions of systems with a = 2.175 and a — 1 as the 
influence of a can be attributed to the effective interac- 
tion c'. 

Next we discuss the momentum distribution for the 
lowest state with K = ±2ir/L. For these cases, unequal 
mass will have significant impact on the wavefunction 
and momentum distribution as the wave functions are not 
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FIG. 4: (Color online) Single particle momentum distribution 
of the ground state for 7 = 5 and a = 1 (dot), 7 = 5 and 
a — 2.175 (solid line), 7 = 500 and a = 1 (square), 7 = 500 
and a — 2.175(dashed line). 



FIG. 6: (Color online) Single particle momentum distribution 
of the lowest state with K = ±2n/L and a — 1 for 7 = 5 and 
7 = 500. 
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momentum distributions become broader and the height 
of the main peak shrinks. As a comparison, in Fig. 6 we 
display the momentum distribution for the lowest states 
of the equal-mass systems with K = ±2tt/L. Due to 
the exchange symmetry of the wave function, the mo- 
mentum distributions ni(p) — n^ip) when mi = 7712. 
For K = 2ir/L, the momentum distribution displays two 
peaks around and 2tt/L, while the distribution has two 
peaks around and —2ir/L for K = —2ti/L. With the 
increase in the strength of the repulsive interaction, the 
height of the peaks also shrinks and the momentum dis- 
tributions become broader. 



V. CONCLUSIONS 



FIG. 5: (Color online) Single particle momentum distribution 
of the lowest state with K = ±2tt/L and a — 2.175 for (a) 
7 = 5 and (b) 7 = 500. 

invariant under the exchange of two heteronuclear atoms. 
Consequently, the momentum distributions ni(p) and 
ri2(p) are different when mi ^ m,2. In Fig. 5, we display 
the momentum distributions for the lowest states of the 
systems with K = ±2n/L and a = 2.175. For K = 2-k/L, 
the momentum shows the peak of m (p) is around and 
that of ri2{p) is around 2it/L. The momentum distribu- 
tion for the noninteracting system with the same K and a 
has a similar structure as the lowest state with K = 2tt/L 
is composed of the product of single atom states charac- 
terized by fci = and k% = 2n/L. For K = —2n/L, 
the momentum distributions show that the peak of n%(p) 
is around and that of U2(p) is around —2ti/L. With 



In summary, we have exactly solved the problem of 
two heteronuclear atoms interacting with a short-range 
8 potential in a ring trap. For the general case with an 
arbitrary mass ration a, the relative momentum is cou- 
pled with center-of-mass momentum due to the periodic 
boundary condition. By taking the generalized Bethe- 
ansatz-type wave function, we have presented a detailed 
derivation of the analytical solution. While a = 1, our 
solution reduces to the Bethe-ansatz solution of the two- 
particle Lieb-Liniger model. We have applied our an- 
alytical results to studying the energy spectrum, den- 
sity distribution and momentum distributions of the het- 
eronuclear system, which differ from the case of identical 
atoms. 

Acknowledgments.— This work is supported by 
the NSF of China under Grants No.10821403 and 
No. 10974234, and National Program for Basic Research 
of MOST. 



6 



[1] S. Inouye, J. Goldwin, M. L. Olsen, C. Ticknor, J. L. 

Bohn, and D. S. Jin, Phys. Rev. Lett. 93, 183201 (2004). 
[2] C. A. Stan, M. W. Zwierlein, C. H. Schunck, S. M. F. 

Raupach, and W. Ketterle, Phys. Rev. Lett. 93, 143001 

(2004). 

[3] C. Ospelkaus, S. Ospelkaus, L. Humbert, P. Ernst, K. 
Sengstock, and K. Bongs, Phys. Rev. Lett. 97, 120402 
(2006). 

[4] F. Ferlaino, C. DErrico, G. Roati, M. Zaccanti, M. In- 
guscio, G. Modugno, and A. Simoni, Phys. Rev. A 73, 
040702R (2006); M. Zaccanti, C. DErrico, F. Ferlaino, 
G. Roati, M. Inguscio, and G. Modugno, Phys. Rev. A 
74, 041605R (2006). 

[5] M. Taglieber, A. C. Voigt, T. Aoki, T. W. Hansen, and 
K. Dieckmann, Phys. Rev. Lett. 100, 010401 (2008). 

[6] E. Wille, et ai, Phys. Rev. Lett. 100, 053201 (2008). 

[7] M. Iskin and C. A. R. Sa de Melo, Phys. Rev. Lett. 97, 
100404 (2006); 

[8] S.-T. Wu, C.-H. Pao, and S.-K. Yip, Phys. Rev. B 74, 

224504 (2006); C.-H. Pao, S.-T. Wu, and S.-K. Yip Phys. 

Rev. A 76, 053621 (2007). 
[9] M. A. Baranov, C. Lobo, and G. V. Shlyapnikov, Phys. 

Rev. A 78, 033620 (2008). 
[10] S. B. Papp and C. E. Wieman, Phys. Rev. Lett. 97, 

180404 (2006). 

[11] A.-C. Voigt, M. Taglieber, L. Costa, T. Aoki, W. Wieser, 
T. W. Hansch, and K. Dieckmann, Phys. Rev. Lett. 102, 



020405 (2009). 

[12] F. Deuretzbacher, K. Plassmeier, D. Pfannkuche, F. 

Werner, C. Ospelkaus, S. Ospelkaus, K. Sengstock, and 

K. Bongs, Phys. Rev. A, 77, 032726 (2008). 
[13] G.-D. Lin, W. Yi, and L.-M. Duan, Phys. Rev. A 74, 

031604(R)(2006). 
[14] D. Blume, Phys. Rev. A 78, 013613 (2008); D. Blume 

and K. M. Daily, Phys. Rev. Lett. 105, 170403 (2010). 
[15] Th. Busch, B. G. Englert, K. Rzazewski, and M. Wilkens, 

Found. Phys. 28, 549 (1998). 
[16] M. A. Cirone, K. Coral, K. Rzazewski, M. Wilkens, J. 

Phys. B, 34, 4571 (2001). 
[17] S. Frankc- Arnold, S. M. Barnett, G. Huyet, and C. Sail- 

liot, Eur. Phys. J. D 22, 373 (2003). 
[18] E. H. Lieb and W. Liniger, Phys. Rev., 130, 1605 (1963). 
[19] P. Shea, B. P. van Zyl and R.K. Bhaduri, Am. J. Phys, 

77, 511 (2009). 

[20] R. T. Piil, N. Nygaard, and K. Molmer, Phys. Rev. A, 

78, 033611 (2008). 

[21] S. Gupta, K. W. Murch, K. L. Moore, T. P. Purdy, and D. 

M. Stamper-Kurn, Phys. Rev. Lett. 95, 143201 (2005); 

A. S. Arnold, C. S. Garvie, and E. Riis, Phys. Rev. A 

73, 041606(R) (2006). 
[22] C. Ryu, M. F. Andersen, P. Clade, V. Natarajan, K. 

Helmerson and W. D. Phillips, Phys. Rev. Lett. 99, 

260401 (2007). 



